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Abstract 

The axial anomaly of Ginsparg- Wilson fermion operator D is dis- 
cussed in general for the operator R which enters the chiral symmetry 
breaking part in the Ginsparg- Wilson relation. The axial anomaly and 
the index of D as well as the exact realization of the Atiyah-Singer 
index theorem on the lattice are determined solely by the topological 
characteristics of the chirally symmetric operator D c in the chiral limit 
R -> 0. 
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Ginsparg- Wilson (GW) relation |I] has been playing an important role in 
recent theoretical developments of the chiral symmetry on the lattice @, [jj. It 
was derived by Ginsparg and Wilson in 1981 as the remnant of chiral symme- 
try on the lattice after blocking a chirally symmetric theory with a chirality 
breaking local renormalization group transformation. The GW relation for 
exactly massless Dirac fermion operator D is 

D l5 + l5 D = 2D l5 RD, (1) 

where the chiral symmetry breaking on the RHS of Eq. ([I]) constitutes an 
invertible hermitian operator R which is local in the position space and trivial 
in the Dirac space. 

In topologically trivial background gauge field, D~ x exists and Eq. (|l]) can 
be written in terms of the fermion propagator D^ 1 

ITS + lsD~ l = 2 l5 R (2) 

This equation displays explicitly that the chirality breaking part in the GW 
fermion propagator is dictated by the local operator R, in contrast to those 
non-local breakings due to adding a chiral symmetry breaking term explicitly 
to the action. 

The GW relation, Eq. (|1|) can be rewritten as 

£>7 5 (1 - RD) + (1 - DR) l5 D = (3) 

Then it becomes evident that the lattice action A = tpDip is invariant under 
the exact chiral transformation on the lattice [|J : 

V> -> exp[^7 5 (l--R-D)]'0, (4) 
i> -> $ exp[iO(l - DR)j B ], (5) 

where 9 is a global parameter. Consequently the axial anomaly can be deduced 
from the change of fermion integration measure under the exact chiral trans- 
formation (f|)-([5p, and its sum over all sites is equal to the index of D which 
is a well defined integer || |J. The Ginsparg- Wilson fermion circumvents the 
Nielson-Ninomiya no-go theorem || by breaking the continuum chiral symme- 
try {75, D} = but keeping the exact chiral symmetry on the lattice, Eq. (|), 
hence D can be constructed to be local and free of species doubling. 

The salient feature of the GW relation Eq. ([I]) in topologically nontrivial 
background gauge field is that its chiral symmetry breaking part does not have 
any effects on the zero modes. This ensures that D has exact zero modes with 
definite chiralities. Consequently the axial anomaly and the index of D are 
invariant with respect to R, and thus they are determined solely by the chirally 
symmetric operator D c in the chiral limit R — > 0, where the GW relation 
is completely turned off. The role of the GW relation is to transform the 
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nonlocal and sometimes singular ( i.e., divergent ) D c into a class of local and 
well defined fermion operators D while keeping the topological characteristics 
of D c invariant. If we have a D c which can possess nonzero indices, then 
any D constructed from this D c using the general solution of GW relation 
has the same topological characteristics. Conversely, if we have a Ginsparg- 
Wilson D which has proper indices in topologically nontrivial gauge fields, 
then there must exist a chirally symmetric D c which has the same topological 
characteristics of D, and we can use this D c to construct a class of GW fermion 
operators. The purpose of this paper is to discuss the i?-invariance of the axial 
anomaly, the index and the Atiyah-Singer index theorem on the lattice, and 
show that their realizations on the lattice are actually beyond the control of 
the GW relation but determined solely by the topological characteristics of 
the chirally symmetric D c in the chiral limit R — > 0. For the axial anomaly 
recently constructed by Liischer [71], we show that the term proportional to FF 
is independent of R but depends on the topological characteristics of D c , while 
the divergence term may depend on R in a very complicated way, but we show 
that its relationship to the divergence of axial vector current is i?-invariant. 

We begin by reviewing the general solution of the Ginsparg- Wilson rela- 
tion which has been discussed in ref. ||. Our purpose here is to set up the 
notations. The general solution of Eq. (|l|) can be formally represented by 

D = D c (l + RDJ- 1 = (l + D.R^D, (6) 

where D c is the chirally symmetric fermion operator of Eq. ([!]) in the chiral 
limit R -»• 0, 

£ C 75 + 75 £ C = 0. (7) 

It is instructive to show that these two different forms of the general solution in 
@ are equal and they satisfy Eq. ([!]). The details are shown in the following. 

D^l + RDc)- 1 
= (1 + D C R)- 1 (1 + D C R)D C (1 + RD C )- 1 
= (1 + D C R)- 1 D C (1 + RD C )(1 + RD C )- 1 
= (l + DcR^Dc 

Using the general solution (|) to substitute D on the RHS of Eq. (0), and 
inserting Eq. ([7J), we obtain 

2D l5 RD 

= 2(1 + D c R)- l D clh RD c {l + RD c y l 

= (l + D c R)~ 1 D c R l5 D c (l + RD c )- 1 + 
(1 + D c R)- l D cl r RD c (l + RD C )- 1 + 
(1 + D c R)-\ l5 D c + D cTo )(l + RDcY 1 
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= (l + D c R)~\l + D c R) l5 D c (l + RD C )- 1 + 

(1 + D c R)~ l D cl5 {l + RD C )(1 + RDc)- 1 
= 75D + D75 

For any D c ( hence D ) satisfying the hermiticity condition 

75 D c75 = Dl (8) 

then D c is antihermitian and it is evident from M that there is one to one cor- 
respondence [] between D c and a unitary operator V satisfying the hermiticity 
condition ( 75^75 = ) 

D e = (l + V)(l - V)-\ V = (D C - 1)(D C + I)" 1 . (9) 

It is interesting to note that V exists for all D c , in particular, when D c — > ±zoo, 
V — > +1; but when D c — > 0, V — > — 1. The eigenvalues of V fall on a unit circle 
with center at zero. On the other hand, D c becomes singular ( with simple 
poles ) when V has real eigenvalues +1. This is exactly what we want it to 
happen when the background gauge field is topologically nontrivial. However 
D is still well defined even when D c is singular. Substituting D c into (||), we 
obtain the general solution of D in terms of a unitary operator V satisfying 
the hermiticity condition [EL 

D = (l + y)[(I- V) + R(l + V)}- 1 (10) 
= [(1-V) + (1 + V)R]-\1 + V). (11) 

The general solution Eqs. (|T0|), (|TT|) can be casted into many different forms, 
for example 0, 

D = ^=(1 + W)—^= (12) 



where W is a unitary operator which can be expressed in terms of the V and 
R. However, since W depends on R, using this form could easily obscure some 
simple facts especially in the chiral limit R — *■ 0. Therefore we refrain from 
using Eq. flT2"| ) for analytic studies involving the general solution of the GW 
relation. 

In general, the chirally symmetric operator D c is constructed such that it 
agrees with the continuum Dirac operator 7 At (<9 At + iA^) in the classical con- 
tinuum limit. It is obvious from Eq. (H) that if D c has species doubling, D 
must have species doubling too. In order to avoid species doubling for a chi- 
rally symmetric operator, D c must be nonlocal. We refer to ref. || for a few 

"'^This observation is made by Herbert Neuberger in a private communication to the 
author. 



3 



explicit examples of D c . Next we consider the topological characteristics of a 
Dirac operator D by comparing its index, index(D) to the topological charge 
Q of the smooth background gauge field. The index of D is defined as iV_ — iV + 
where N + (N_) denotes the number of zero modes of D with positive ( nega- 
tive ) chirality. If index (£)) = for any smooth gauge configurations, then D 
is called topologically trivial. If index(D) = Q, then D is called topologically 
proper. If D is not topologically trivial but index(D) ^ Q, then D is called 
topologically improper. Examples of these three different cases are given below 



and detailed discussions in ref. |]10| . From Eq. (y), it is obvious that any 
zeromode of D is a zeromode of D c and vice versa. Therefore the index of D 
is the same as the index of D c , independent of R. 

index[L>, A] = index[£> c , A] = N_[D C , A] - N + [D C , A] (13) 

where we have indicated explicitly that the number of zero modes and the 
index are functionals of D c ( D ) and the gauge field A. However, if D c is well 
defined, its index must be zero. This can be easily proved in the following. 
Proof : Consider R = rl, then index(D) = r J2 X ^ r [l5D(x, x )} S HI- m the 
limit r — > 0, D — > D c , if D c is well defined, then J2x tr^D^x, x)] is finite, thus 
multiplying r gives index(D) = 0. This completes the proof. 
Therefore in order for D c to possess non-zero indices, D c must be singular 
( or equivalently, V has real eigenvalue pairs ±1 ) in topologically nontrivial 
background gauge field || . 

To summarize above discussions, we list the necessary requirements for D c to 
enter Eq. (^) for constructing a local and well defined D such that D could 
possibly reproduce the continuum physics : 

(a) D c is chirally symmetric and agrees with 7^(9^ + iA^) in the classical 

continuum limit. 

(b) D c is free of species doubling. 

(c) D c is nonlocal. 

(d) D c is well defined in topologically trivial background gauge field. 

(e) D c is singular ( or equivalently, V has real eigenvalue pairs ±1 ) in topo- 

logically non-trivial background gauge field. 

It should be noted that each one of these constraints is required to be satisfied. 
For example, if (b) is not satisfied, then the index of D c must be incorrect no 
matter other constraints are satisfied or not. The most difficult task is to find 
a D c satisfying (e). 

At this point, it is instructive to consider examples of topologically proper, 
trivial and improper D c respectively. First we consider the Neuberger-Dirac 
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operator [IT] in the chiral limit ( R — > 



Dc = 2M f3V' V^A^A,)" 1 ' 2 (14) 

where M is an arbitrary mass scale, D w is the Wilson-Dirac fermion operator 
with negative mass —mo and Wilson parameter r w > 

A, = -^o + ^[7m(V* + V„) -r w V*V M ] (15) 

where V M and V* are the forward and backward difference operators defined 
in the following, 

V^{x) = U IJi {x)'4){x + ft) - ip(x) 
V*i/j(x) = ip(x) - Ul(x - fi)ip(x - fi) 

Using Eq. (f|), we obtain the generalized Neuberger-Dirac operator 

D = 2M(1 + V)[(1-V) + 2MR(1 + V)}- 1 (16) 

For m G (0, 2r w ), both D and .D c are topologically proper and reproduce the 
exact index theorem on a finite lattice for smooth background gauge fields as 



first demonstrated in 13, 14 . On the other hand, for mo < 0, that is, D w with 



a positive mass, then D and D c becomes topologically trivial ITI3]. We note 



that in both trivial and proper cases D is free of species doubling and is local 
for properly chosen R. For m > 2r w , both D and D c are either topologically 
improper or trivial. We emphasize that these ranges of mo values depend on 
the gauge configurations, and the ranges given above are exact only in the 
free field limit. We refer to ref. [JD| for further discussions on the topological 
characteristics of the Neuberger-Dirac operator. 

It is interesting to note that the GW relation does not play any role in 
determining the index of D c . The GW relation only ensures that the chiral 
symmetry is broken in such a delicate way that the exact zero modes of D c 
are also exact zero modes of D with the same chiralities, and thus the index 
of D is equal to the index of D c . The role of R is to convert the nonlocal and 
sometimes singular D c into a local and well defined D, as well as to determine 
the locality of D. If one could not find any D c satisfying all constraints (a)-(e), 
then the GW relation does not have much use at all. On the other hand, if one 
knows one D c satisfying all constraints (a)-(e), then the resulting D must have 
exact chiral symmetry on the lattice as well as all attractive features pointed 



out in 15, 16 



The fermion integration measure in general is not invariant under the exact 
chiral transformation (||)-([5j), and this gives the axial anomaly on the lattice 

i 

q(x,A;D)=tr[ l5 (RD)(x,x)} (17) 
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where tr denotes the trace over Dirac indices. We note that q(x, A; D) is a 
function of x and the gauge field in the neighborhood around x, and a 
functional of D. The sum of q(x, A; D) over all lattice sites is shown [|, |TT[ 
to equal to the index of D 

J2q(x,A;D) = £tr [75 x)] 

X X 

= AL [D, A] - N + [D, A] = indexLD, A] (18) 

where N + (N_) denotes the number of zero modes of D with positive ( negative 
) chirality. In this paper, we restrict our discussions to U(l) gauge theory with 



single fermion flavor. Since the index of D is independent of R [ Eq. (13) ], 
the total axial anomaly must be invariant under any smooth deformations of 
R 

5 R J2q(x,A;D) = (19) 

X 

Recently Luscher |7j has proved that for R = 1/2, if the axial anomaly 
q(x, A; D) satisfies 

J2Wx,A;D)=0 (20) 

X 

for any local deformations of the gauge field, then 

q(x, A; D) = 7 c[D] e^ vX oF^{x)F Xa (x + jl + u) + d^G^x, A; D) (21) 

where 7 is a constant, c[D] is a functional of D which only depends on the 
topological characteristics of D, and 

dpG^x, A;D) = J2 [G„(x, A; D) - G^(x - ft, A; D)] . (22) 

The explicit form of the current G^(x,A;D) is supposed to be very compli- 
cated. We note that c[D] is not shown explicitly in 0. Since the total axial 
anomaly is invariant under smooth deformations of R, Eq. (|T9"D, as well as for 
any local deformations of the gauge field, Eq. (pO"D , then the theorem proved 
by Luscher 0] can be applied to the case of the general R 0, and the axial 
anomaly can be written in the same form as Eq. ( pip on an infinite lattice. 
But on a finite lattice, q(x, A; D) must be modified by the boundary effects. 
However, if D is local, the boundary effects are negligible, then Eq. (pip can 
be used for the axial anomaly on a finite lattice. Since the locality of D is 
controlled by R, we can infer that all boundary effects must have dependence 
on R and hence they do not enter the FF term. [ As we will show later, the 
FF term does not depend on R. ] Furthermore they must be in the form of 

2 This is confirmed by Martin Luscher in a private communication to the author. 
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the divergence term, otherwise Eq. fl21| ) cannot be consistent with Eq. (f2§|). 
the divergence of axial vector current on a finite lattice. Therefore we conclude 
that Eq. Q2"T| ) is applicable for a finite lattice as well as for an infinite lattice, 
and the boundary effects only enter the divergence term. Now we proceed 
our discussions for a finite lattice embedded in a 4-dimensional torus. After 
summing q(x, A; D) over all sites, the divergence term is removed, then Eq. 
(PD gives 

7 c[D) tnv\aF^{x)F Xa (x + p. + i>) 

x 

= N_[D, A] — N + [D,A] — index[D, A] (23) 



Since index[D, A] is independent of R [ Eq. (|T^) ], then Eq. fl23|) implies that 
c[D] must be also independent of R, 

c[D] = c[D c ] (24) 

So we have proved that all R dependence in q(x, A; D) only resides in the 
current G ^(x , A; D) . Then Eq. (p3|) can be rewritten in the following. 

7 c[D c ] ^ ^^\ctF^ v (x)Fx ct {x + fx + v) 

X 

= N_[D C , A] - N + [D C , A] = mdex[D c , A] (25) 

Since index[D c , A] = if D c is topologically trivial, then c[D c ] also vanishes for 
trivial D c . Now we consider two topologically proper D c which are identical 
except their electric charges of the fermion field. Then for a given topologically 
non-trivial gauge configuration, Eq. ( p5| ) implies that the ratio of c[D c ]'s corre- 
sponding to these two D c must equal to the ratio of two integers corresponding 
to their indices. It follows from this that c[D c ] must be an integer. Therefore 
we can fix c[D c ] = 1 for electric charge equal to one and then determine the 
value of the constant 7 by evaluating J2x e fj.uXaF fJiU (x)F\ a (x + (l + v) on the 
LHS for the simplest nontrivial gauge configuration with constant field ten- 
sor, namely, F 12 = 27i/(LiL 2 ) and F 34 = 27r/(L 3 L4) while other F's are zero. 
After restoring the unit of electric charge, we obtain 7 to be integer multiple 
of e 2 /(327r 2 ) for U(l) gauge theory with single fermion flavor. For Neuberger- 
Dirac fermion operator [11]], this integer is determined to be unity by a recent 



perturbation calculation of q(x,A;D) JT^|. For the general Ginsparg- Wilson 



fermion with R = rl, and D c satisfying constraints (a)-(e) Q, our recent per- 
turbation calculation ]Tj| of q(x, A; D) show that 7 is independent of r and 
equal to e 2 /(327r 2 ). Then Eq. ( |25| ) becomes 

e 2 
32tt 2 



c[D e ] J2 ^uX*F^(x)F Xa (x + fi + u) = N4D C , A] - N + [D C , A] (26) 



3 In rcf. |L8), it is only required that the chiral fermion operator D c (p) is free of species 
doubling and in the free field limit behaves like ij^p^ as p — > 0. 
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where 



1 if D c is topologically proper 

c[D c ] — { if D c is topologically trivial (27) 

integer ^ 0, 1 if D c is topologically improper 



For any chirally symmetric D c satisfying constraints (a)-(e), Eq. ^2q ) is the 
exact realization of the Atiyah-Singer index theorem on the lattice and it holds 
for any Ginsparg- Wilson fermion operator D = D c (l+ RD^ 1 . 

Although the explicit form of the current G M (x, A; D) in Eq. ( pT|) is sup- 
posed to be very complicated, however, it is instructive to compare Eq. fl2T|) 
to the divergence of axial vector current on the lattice III H, 



AT, N- 



d,(4(x, A- D)) = 2 q(x, A; D) + 2 YM( X W+( X ) - 2j2[cf>t(x)H t _(x) (28) 

s t 

where <j) s + and are normalized eigenfunctions of D ( D c ) with eigenvalues 
X s = X t = and chiralities +1 and —1 respectively. Eliminating q(x, A; D) 
from both equations, we obtain 

1 e 2 
d M G M (x, A; £>) = -d^Jfe, A-D))- — c[D] e^ Xa F^{x)F Xa (x + fi + v) 

N+ N- 

-E^+^lV+C^ + E^-C^lV-^) (29) 

s i 

It is interesting to note that all dependences on R are in the currents G M and 
J^. The zero modes and the FF term do not depend on R as we have discussed 
above. Then we obtain 



G,{x,A-D)- l -(Jl{x,A-D)) 



(30) 



This equation provides a ^-invariant constraint between the current 
G^x, A; D) and the axial vector current J^(x, A; D). The axial vector current 
J^(x, A; D) is defined |l[ in terms of the kernel K^(x,y, z, A; D) through the 
equation 

4(x,A;D) =Y / ^+yKl(x,y,z,A; D)^ x+Z (31) 
where K\ is expressed in terms of the kernel for the vector current, 



K= o (^75 -75^)- (32) 



Ginsparg and Wilson [Q proved that the kernel for the vector current 

K^x, y, z, A; D) is equal to D(x + y, x + z, A) times the sign of (y — z)^ and 
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times the fraction of the shortest length paths from x + z to x + y which pass 
through the link from x — jl to x. Finally the axial anomaly in Eq. fl2l|) 
becomes 

e 2 

q(x, A; D) = — — c[D c ] e lluX aF fM/ (x)F Xa (x + p, + v) + <9 M G M (x, A; D) (33) 
62n 

where c[D c ] is defined in Eq. ( f27f) and the divergence term d^G^(x, A; D) is 



related to the divergence of axial vector current by Eqs. ( |29| ) and (|30"D. 

To summarize, we have clarified the role of the Ginsparg- Wilson relation in 
the realization of the axial anomaly and the exact Atiyah-Singer index theorem 
on the lattice. The crucial point is the existence of a topologically proper 
chirally symmetric Dirac operator D c in the chiral limit where the Ginsparg- 
Wilson relation is turned off. If D c satisfies the constraints (a)-(e) discussed in 
this paper, then any Ginsparg- Wilson fermion operator D constructed by Eq. 
(^) could reproduce the correct continuum physics, in particular, the correct 
axial anomaly and the exact index theorem. From this point of view, the 
underlying reason why Neuberger-Dirac fermion really works is because the 
existence of the topologically proper D c ( or V ) [ Eq. (|14|) ] which stems from 
the Overlap formalism |T9|, The Ginsparg- Wilson relation only provides a 
scheme to break the continuum chiral symmetry in such a delicate way that 
the exact chiral symmetry on the lattice is preserved and the axial anomaly 
and the index are invariant while smooth deformations of R away from zero 
bring the non-local and sometimes singular D c into a class of local and well 
defined Ginsparg- Wilson fermion operators D. The most crucial and highly 
non-trivial task is to find a 'good seed' D c which at least satisfies all constraints 
(a)-(e). So far there is only one good D c found explicitly. However, there is 
no indications that D c is unique and thus the search for computationally less 
expensive D c is highly desirable. Since the Overlap formalism did provide 
vital insights in discovering a good D c JTl||, it would be interesting to see 



whether it could guide us to find a better one. Mathematically, it is interesting 
to understand what actually determines the topological characteristics of an 
operator and how to construct a topologically proper operator in a systematic 
way. For the axial anomaly ([H]) recently constructed by Luscher J7|, we show 
that the FF term is i?-invariant and the proportional constant is e 2 /(327r 2 ) 
times c[D c ] which only depends on the topological characteristics of D c \ and 
the divergence term is related to the divergence of axial vector current by an 
/^-invariant constraint, Eqs. (E^) and For any chirally symmetric D c 

satisfying constraints (a)-(e), the Atiyah-Singer index theorem can be realized 
exactly on the lattice, Eq. Q2BD, and it holds for any Ginsparg- Wilson fermion 
operator D = D c (l + FLD c )~ l . 
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